Exact solution for mean first-passage time on a pseudofractal scale-free web 
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The exphcit determinations of the mean first-passage time (MFPT) for trapping problem are 
limited to some simple structure, e.g., regular lattices and regular geometrical fractals, and deter- 
mining MFPT for random walks on other media, especially complex real networks, is a theoretical 
challenge. In this paper, we investigate a simple random walk on the the pseudofractal scale-free 
web (PSFW) with a perfect trap located at a node with the highest degree, which simultaneously 
exhibits the remarkable scale-free and small-world properties observed in real networks. We obtain 
the exact solution for the MFPT that is calculated through the recurrence relations derived from 
the structure of PSFW. The rigorous solution exhibits that the MFPT approximately increases as a 
power-law function of the number of nodes, with the exponent less than 1. We confirm the closed- 
form solution by direct numerical calculations. We show that the structure of PSFW can improve 
the efficiency of transport by diffusion, compared with some other structure, such as regular lat- 
tices, Sierpinski fractals, and T-graph. The analytical method can be applied to other deterministic 
networks, making the accurate computation of MFPT possible. 

PACS numbers: 05.40.Fb, 89. 75. He, 05.60.Cd, 89. 75. Da 



I. INTRODUCTION 

The problem of random walks (diffusion) is a major 
theme and understanding its behavior is central to a 
wide range of applications P, 0, B 0| • Among various 
interesting issues of random walks, trapping plays a sig- 
nificant role in an increasing number of disciplines, e.g., 
physics d, 0] , society , computer |8i] , and biology Q , 
to name a few. The classic trapping issue first intro- 
duced in [lo| is a random- walk problem, in which a trap 
is placed at a given location, absorbing all walkers that 
visit it. 

An important quantity related to trapping problem is 
the trapping time (first-passage time, survival time, or 
the mean walk length) or the mean time to absorption. 
The trapping time of a given site s is the expected time 
for a walker starting from s to first reach the trap. This 
quantity is useful in the stud y o f transport-limited reac- 
tions [nl, [l^l , target search [13l , [T^ and other physical 
problems. The average trapping time, also known as the 
mean first-passage time (MFPT), characterizes the pro- 
cess of trapping. It is defined as the average of survival 
times over all starting sites. 

In the past few decades, there has been considerable in- 
terest in computing the mean first-passage time, in order 
to obtain the dependence of this primary quantity on the 
system size or other parameters. In a seminal work, us- 
ing the approach based on generating functions, Montroll 
derived the rigorous results for MFPT of random walks 
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on regular lattices with a variety of dimensions [l^ . Re- 
cently, by applying a decimation procedure, Kozak and 
Balakrishnan obtained the accurate solutions for MFPT 
on a family of Sierpinski fractals [Tsl . [l^ ; using an analo- 
gous but different method, Agliari got the exact expres- 
sion for the MFPT for a random walker on T- fractal [l3| ■ 
In spite of these rigorous results, the explicit determina- 
tion of MFPT for random walks with a trap on other 
media is still open [isf. 

It is well established that the scaling of MFPT is 
related to the underlying structural properties of the 
media in which the walkers are confined [l|, 0, d, 
Extensive empirical studies [l^, have revealed that 
most real networked systems share some striking fea- 
tures, such as scale-free behavior [2T] and small-world 
effects ^2^. These newly- found properties have a pro- 
found effect on almost all dynamical processes taking 
place on the networks [23l. HI, fsB'l, including disease 
spreading games [23, synchronization and so 
on. Very recently, a lot of activities have been devoted 
to the study of influences of scale-free and small-world 
characteristics on the behavior of random walks, uncov- 
ering many unusual and exo tic p henomena about random 
walks [2i,l3i[3l|,|33,[3l[33,[3l[3i,[33. However, to best 
of our knowledge, rigorous solution for MFPT of trapping 
problem on scale-free small- world networks is missing. 

In this paper, we study the classic trapping problem on 
a deterministic network, called pseudofractal scale-free 
web (PSFW) [111. We focus on a peculiar case with the 
trap fixed at a hub node (node with the highest degree). 
The PSFW is a very useful toy model that captures si- 
multaneously scale-free small- world properties, thus pro- 
vides a good facility to investigate analytically trapping 
process upon it. We derive an exact formula for the 
mean first-passage time characterizing the trapping pro- 
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From Eq. ([2]), one can easily see that at each step the 
degree of a node doubles, i.e., 



FIG. 1: The first three generations of the pseudofractal scale- 
free web. 



cess. The analytic approach is based on an algebraic it- 
erative procedure obtained from the particular construc- 
tion of PSFW. The obtained rigorous result shows that 
the MFPT grows as a power-law function of the num- 
ber of network nodes with the exponent less than 1, 
which implies that in contrast to regular lattices, Sier- 
pinski fractals, and T-graph, the PSFW tends to speed 
up the diffusion process. Our study opens the way to 
theoretically investigate the MFPT of a random walker 
on a wide range of deterministic networks [s^ liO, lU, \^ . 



II. THE PSEUDOFRACTAL SCALE-FREE WEB 

Here we introduce the pseudofractal scale-free web de- 
fined in a recursive way |38| . which has attracted an 
amount of attention [H, |44|. l45l l46j . We investigate the 
PSFW model because of its intrinsic interest and its de- 
terministic construction, which allows one to study ana- 
lytically their topological properties and dynamical pro- 
cesses on it. 

The pseudofractal scale-free web, denoted by Gn after 
n (n > 0) generation evolution, is constructed as follows. 
For n = 0, Go is a triangle of edges connecting three 
nodes (vertices, sites). For n > 1, Gn is obtained from 
Gn-i- every existing edge in Gn-i introduces a new node 
connected to both ends of the edge. Figure [1] illustrates 
the construction process for the first three generations. 

According to the network construction, one can see 
that at each step Ui (rii > 1) the number of newly in- 
troduced nodes is L{ni) = 3"*. From this result, we can 
easily compute the network order (i.e., the total number 
of nodes) Vn at step n: 



Vn 



(1) 



ni=0 



Let ki{n) be the degree of a node i at time n, which 
entered the network at step rii [rii > 0). Then 



ki(n) — 2ki{n — 1). 



(3) 



(2) 



The PSFW presents some typical characteristics of 
real-life networks in nature and society. It has a power- 
law degree distribution P(fc) ^ k^'^ with the exponent 
7 = 1 + 1^ [38]. Its average path length (APL), 
defined as the mean of shortest distance between all 
pairs of nodes, increases logarithmically with network or- 
der [H, [4^. In the large network order limit, the aver- 
age clustering coefficient tends to |. Thus, the PSFW 
exhibits small- world behavior [2^ . In addition, its node 
degree correlations are negative, the average degree of 
nearest neighbors, fc„„(A:), for nodes with degree k is ap- 
proximately a power-law function of k with a negative 
exponent -(2- j^^) [44]. 

While the graph introduced above visually looks very 
similar to a fractal (see Fig. [1]), the similarities are only 
superficial. The distinction is obvious. Fractals have 
a finite dimension ( i.e., their APL grows as a power of 
the vertex number) 47]. In sharp contrast, the network 
under consideration has an infinite dimension (i.e., its 
APL increases slower than any power of the network or- 
der) [48j . This is why it is called pseudofractal scale-free 
web. 

After introducing the PSFW, in what follows we will 
study the mean first-passage time for random walks with 
a single immobile trap on the web, by applying a method 
similar to but different from that introduced in llSl lid] . 



III. FORMULATION OF THE PROBLEM 

In this section we formulate the problem of a simple 
unbiased Markovian random walk of a particle on PSFW 
Gn in the presence of a trap or a perfect absorber located 
on a given node. In order to distinguish different nodes, 
we label all the nodes belonging to G„ in the following 
way. The initial three nodes in Gg are labeled as 1, 2, 
and 3, respectively. In each new generation, only the 
new nodes created at this generation are labeled, while 
the labels of all old nodes remain unchanged, i.e., we 
label new nodes as M + 1, M + 2, . . . , M + AM, where 
M is the total number of the pre-existing nodes and AM 
is the number of newly-created nodes. Eventually, every 
node is labeled by a unique integer, at time n all nodes 
are labeled from I to Vn — - — 1 see Fig. [TJ 

We locate the trap at node 1 [i^, denoted as it- Note 
that the particular selection we made for the trap loca- 
tion makes the analytical computation process (that will 
be shown in detail in the next section) easily iterated as 
we can identify the trap node ir since the first genera- 
tion. At each time step (taken to be unity), the particle, 
starting from any node except the trap zt, moves to any 
of its nearest neighbors with equal probability. It is eas- 
ily seen that in the presence of the trap ix fixed on node 
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1, the walker (particle) will be inevitably absorbed [32| . 
This random walk process can be described by a specify- 
ing the set of transition probabilities Wij for the particle 
of going from node i to node j. In fact, all the entries 
constitute a matrix W that is a {Vn — l)-order sub-matrix 
of (Z"^)A with the first row and column corresponding 
to trap being removed, where Z and A are separately the 
adjacency matrix and diagonal degree matrix of G„ [50| . 

Let Ti be survival time for a walker initially placed 
at node i to first reach the trap ir- Then, the set of 
this interesting quantity obeys the following recurrence 
equation [5l| 

T, = ^M^,,T, + 1, (4) 
j 

where i ^ ir- Eq- (21) expresses the Markovian property 
of the random walk, it may be recast in matrix notation 
as 

T = WT + e, (5) 

where T = {T2,T3, ■ ■ • ,ry„)^ (the superscript T of the 
vector represents transpose) is a {Vn — l)-dimensional 
vector, e is the {Vn — l)-dimensional unit vector 
,1)^, and W is the transition matrix. From 
Eq. (O, we can easily obtain 

T = Le, (6) 

where 

L={I-W)-\ (7) 

in which I is an identity matrix with an order {Vn — 
1) X {Vn — 1). Equation ([7]) is the fundamental matrix 
of the Markov chain representing the unbiased random 
walk. Actually, the matrix I — W in Eq. ^ is the nor- 
malized discrete Laplacian matrix of G„ whose first row 
and column that correspond to the trap node have been 
deleted. 

Then, the mean first-passage time (MFPT), or the av- 
erage of the mean time to absorption, (T)„, which is the 
average of Ti over all initial nodes distributed uniformly 
over nodes in Gn other than the trap, is given by 

= E = ^ E i: (8) 

•'n ^ „ ^ „ „ 

1—2 1—2 J— 2 

Equation ^ can be easily explained from the Markov 
chain representing the random walk. In fact, the entry 
Lij of the fundamental matrix L for the Markov process 
denotes the expected number of times that the process 
is in the transient state j, being started in the transient 
state i. 

The quantity of MFPT is very important since it mea- 
sures the efficiency of the trapping process: the smaller 
the MFPT, the higher the efficiency, and vice versa. 
Equation ([5]) shows that the problem of calculating (T)„ 



TABLE I: The average of the mean time to absorption ob- 
tained by direct calculation from Eq. Since for large 
networks, the computation of the MFPT from Eq. ([8} is 
prohibitively time and memory consuming, we calculate the 
MFPT for only the first several generations. 
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19699/365 


6 


1095 


117401/1094 
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3282 


702091/3281 


8 


9843 


4205729/9842 



is reduced to finding the sum of all elements of matrix 
L. Notice that the order of L is (F„ — 1) x {Vn — 1), 
where Vn increases exponentially with n, as shown in 
Eq. ID). So, for large n, it becomes difficult to obtain 
(T)„ through direct calculation from Eq. ([5|), one can 
compute directly the MFPT only for the first several gen- 
erations, see TablelD However, the recursive construction 
of PSFW allows one to compute analytically MFPT to 
achieve a closed-form solution, the derivation details of 
which will be given in next section. 



IV. EXACT SOLUTION FOR MEAN 
FIRST-PASSAGE TIME 

Before deriving the general formula for MFPT, (r)„, 
we first establish the scaling relation dominating the evo- 
lution of T" with generation n, where T" is the trapping 
time for a walk originating at node i on the nth genera- 
tion of PSFW. 



A. Evolution scaling for trapping time 

We begin by recording the numerical values of T". Ob- 
viously, for all n > 0, T" = 0; for n = 0, it is a trivial case, 
we have T§ = T^ = 2. For n > 1, the values of can be 
obtained straightforwardly via Eq. ([H]). In the generation 
n = 1, by symmetry we have = Tj^ = 4, = = 3, 
and Tg = 5. Analogously, for n = 2, the solutions are 
Ti = Ti = 8, T| = Tl ^ 6, r| = 10, T? = r| = 4, 
Ti = Tf^ = 8, = Ti, = 10, Ti^ = Ti^ = 5, and 
Ti^ = 9. Table M lists the numerical values of Tp for 
some nodes up to n = 6. 

The numerical values listed in Table |TT] show that for 
a given node i we have T^+i = 2T". That is to say, 
upon growth of PSFW from n to generation n + 1, the 
mean time to first reach the trap doubles. For example, 

T| = 2r| = 4T2' = 8r| = i6T| = 32T2I = 64rO = 
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TABLE II: Mean time to absorption T" for a random walker 
starting from node i on PSFW for various n. Notice that 
owing to the obvious symmetry, nodes in a parenthesis are 
equivalent, since they have the same trapping time. All the 
values are calculated straightforwardly from Eq. ((8}. 
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(7,8) (9,12) (10,11) 


(13,14) (15) 
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160 
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160 


80 144 



128, T| = 2r| = 4T| = 8T| = 16 T| = 32 = 96, 
= 2T| = AT^ ^8T^ ^ 16T^ ^ 64, and so on. This 
is a basic character of random walks on the PSFW, which 
can be established from the arguments below. 

Consider an arbitrary node i in the PSFW G„ after n 
generation evolution of the network. From Eq. ([2|), we 
know that upon growth of PSFW to generation n + 1, the 
degree fc^ of node i doubles. Let the mean transmit time 
for going from node i to any of its hi old neighbors be 
Y \ and let the mean transmit time for going from any of 
its ki new neighbors to one of the ki old neighbors be Z. 
Then we can establish the following underlying backward 
equations 



Y = i + i(l + Z), 

z = | + |(i + y), 



(9) 



which leads to Y = 2. That is to say, the passage time 
from any node i (i G Gn) to any node j {j £ G„) increases 
by a factor of 2, upon the network growth from generation 
n to generation n + 1. Thus, we have T\"+i — 2 T", which 
will be useful for deriving the formula for the mean first- 
passage time in the following text. 



B. Formula for the mean first-passage time 

Having obtained the scaling of mean transmit time for 
old nodes, we now determine the average of the mean 
time to absorption, aiming to derive an exact solution. 
We represent the set of nodes in Gn as A„ , and denote the 
set of nodes created at generation rt by A„. Thus we have 
An = A„ U A„_i. For the convenience of computation, 
we define the following quantities for m < n: 



rpn 

m, total 



E 



(10) 



Then, we have 



total 



rj-rn 



— 1. total 



T- 



.total ■ 



(12) 



Next we will explicitly determine the quantity T!^ total- 
To this end, we should firstly determine r„ tot^i- 

We examine the mean time to absorption for the first 
several generations of PSFW. In the case of n = 1, by 
construction of the PSFW, it follows that = 5(1 + 



n 



\{l + Tl) + \{l+Tl),&ndT^^ 



n) 



Thus, 



1, total 



ieKi 



3 + {Tl + Tl 



n) 



3 + T, 



0. total 



(13) 



Similarly, for n ^ 2 case, = 5(1 + Tf) + i(l + r|). 



rp2 _ 
^8 - 

^ 10 



1 + T?)- 



-Tl) 

I) 



-n 



^9 



Tl), so that 



^(1 + 7^1)- 

T|),r2, = i(i+r|)- 



-Ti) 



i(l + T|),andTf5 = i(l+T|) + i(l- 



2, total 



15 

E^' 

i=7 



= 3' 
= 32 



2 (Tf 
-2 



r| + T|) + (r| + r| + r|) 



2t; 



0, total ■ 



T 



1, total ' 



(14) 



Proceeding analogously, it is not difficult to derive that 

+ ^n-l,total "I" 2 T„_2 total + ■ ■ ■ 
+ 2" ^ total + 2" ^ Tq. total; (15) 



^ n, total — 'J 



and 



T, 



i+i 

n+l, total 



?™+l 



T, 



i+l 



total 
2n-l 



2T 



^ 1, total 



n—1, total . . ■ 

2"r"total, (16) 



where 3" and 3"^^ are indeed the numbers of nodes gen- 
erated at generations n and n + l, respectively. Equa- 
tion minus Eq. times 4 and making use of the 
relation 7^"+^ = 2Tp, one gets 



7p"+l _ on+1 

n+l, total 



t: 



i-i-i 



4 t: 



n, total 

which may be rewritten as 

6 ^n, total 



.total 



1. total 



3"), (17) 



(18) 



and 



total 



= E ^" 



ieAr, 



(11) 



Using T^ j.^^^^ = 11, Eq. (fTSj) is solved inductively 

—n 5 



n I on—1 



n, total — X + 



(19) 



5 



Substituting Eq. ^ for T" ^^^^j into Eq. we have 



n, total n — 1, total ' ^ 



+ _ X 6" + T-^ 
5 



(20) 



Considering the initial condition ^^^^j = 4, Eq. ((20l) is 
resolved by induction to yield 



fn'.total — - X 6" + 3" 



-)n — 1 



(21) 



Plugging the last expression into Eq. ([8]), we arrive at 
the accurate formula for the average of the mean time to 
absorption at the trap located at node 1 on the nth of 
the pseudofractal scale-free web: 



{T)r 



1 



1 



t—Z 

5 X 6" + 2 X 3" + 2"^ 



rpn 



total 



3«+i _^ I 



(22) 



We have checked our analytic formula against numerical 
values quoted in Table |T1 For the range of < rt < 
8, the values obtained from Eq. ((22l) completely agree 
with those numerical results on the basis of the direct 
calculation through Eq. ([8]) . This agreement serves as an 
independent test of our theoretical formula. 

We continue to show how to represent MFPT as a func- 
tion of network order, with the aim of obtaining the scal- 
ing between these two quantities. Recalling Eq. ([T|), we 
have 3"+i = 2Vn - 3 and n + 1 = log3(2K - 3). Hence, 
Eq. ([22]) can be recast as 



i(2K 



{T)r 



For large network, i.e., Vn 



3)i+M + |(2K-3) + M2K-3)i^« 



2Vn 



(23) 



oo. 



(T}„~(K)^, (24) 

where the exponent ||-| < 1. Thus, in the large limit of 
network order Vn, the MFPT increases algebraically with 
increasing order of the network. 

The above scaling of the MFPT with network order 
is different from those previously obtained for other me- 
dia. For example, on regular lattices with large order 
iV, the leading behavior of MFPT (T) is (T) - 
(T) - A^lniV, and (T) ~ N for dimensions d=l,d = 2, 
and d = 3, respectivel y [lol . Again for instance, onpla- 
nar Sierpinski gasket [15(1 ^^'^ Sierpinski tower [IB] in 
3 Euclidean dimensions, the asymptotic behavior scales 
separately as (T) - N^-'^^^ and (T) - N^-^^^, At last, 
on T-fractal, the leading asymptotic scaling behaves as 
(T) - N^ esi ^ Thus, in contrast to regular lattices, 
Sierpinski fractals, and T-fractal, the trapping process on 
the pseudofractal scale-free web is more efficient. It is ex- 
pected that the efficiency of trapping process on stochas- 
tic scale- free networks is also high, since they have similar 
structural features as the PSFW. 



Why is the MFPT for the PFSW far smaller than that 
for other lattices? We speculate that the heterogeneity 
of the pseudofractal scale-free web is responsible for this 
distinction, which may be seen from the following heuris- 
tic argument. In the PFSW, there are a few nodes with 
large degree that are connected to most nodes in the 
web, which results in the logarithmic scaling of the av- 
erage path length with network order (38 . ,46.] . A walker 
starting from some node will arrive at 'large' nodes with 
ease. Since 'large' nodes, including the trap node, are 
linked to one another, so the walker can find the trap 
in a short time. While for other regular lattices, they 
are almost homogeneous, and their average path lengths 
are much larger than that of the PFSW. This leads to a 
long MFPT. It should be noted that we only give a pos- 
sible explanation for the shorter MFPT on the PFSW, 
the genuine reason for this deserves further study in the 
future. 



V. CONCLUSIONS 

In summary, we have investigated the classic trapping 
problem on a deterministically growing network, named 
pseudofractal scale-free web (PSFW) that can reproduce 
some remarkable properties of various real-life networks, 
such as scale-free feature and small-world behavior, and 
thus can mimic some real systems to some extent (to 
what extent it does is still an open question). With the 
help of recursion relations derived from the structure of 
PSFW, we have obtained the solution for the MFPT for 
random walks on PSFW, with a trap fixed at a hub node. 
The exact result shows that the MFPT increases as a 
power-law function of network order with the exponent 
less than 1, which is in contrast to the well-known previ- 
ously obtained results that for regular lattices, Sierpinski 
fractals, and T-graph with order , their MFPT (T) be- 
haves as (T) ~ A^" with a > 1. Therefore, the structure 
of PSFW has a profound impact on the trapping problem 
on it. To the best of our knowledge, our result may be 
the first exact scaling about MFPT for random walks on 
scale-free small- world networks. 

We should stress that although we have only computed 
the MFPT for a particular deterministic network, our an- 
alytical technique could guide and shed light on related 
studies for other deterministic networks by providing a 
paradigm for calculating the MFPT. Moreover, since ex- 
act solutions can serve for a guide to and a test of approx- 
imate solutions or numerical simulations, we also believe 
our accurate closed-form solutions can prompt related 
studies on stochastic networks. At last, as future work, 
it is interesting to compute higher moments of trapping 
time for the PSFW and compare the scaling with that 
of homogeneous fractal lattices [s^. Another future job 
of interest is to study the case when the trap is mobile 
instead of being fixed at one of the hub nodes of the web. 
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